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Abstract. The LHCf experiment has taken data in 2009 and 2010 p-p colhsions at LHC at = 0.9 TeV 
and ^|s - 1 TeV. The measurement of the forward neutral particle spectra produced in proton-proton collisions 
at LHC up to an energy of 14 TeV in the center of mass system are of fundamental importance to calibrate the 
Monte Carlo models widely used in the high energy cosmic ray (HECR) field, up to an equivalent laboratory 
energy of the order of 10^^ eV. 

In this paper the first results on the inclusive photon spectrum measured by LHCf is reported. Comparison of this 
spectrum with the model expectations show significant discrepancies, mainly in the high energy region. 
In addition, perspectives for future analyses as well as the program for the next data taking period, in particular 
the possibility to take data in p-Pb collisions, will be discussed. 



1 Introduction 

The LHCf experiment at LHC has been designed to cali- 
brate the hadron interaction models used in High Energy 
Cosmic Ray (HECR) Physics through the measurement of 
the forward neutral particle produced in p-p interactions. 
The uncertainty caused by the poor knowledge of the in- 
teraction between very high energy primary cosmic rays 
and the earth's atmosphere prevents the precise deduction 
of astrophysical parameters from the observational data, 
thus causing the most important source of systematic error 
in HECR dedicated experiments. Dedicated extensive air 
shower experiments are in place since many years and have 
strongly contributed to our understanding of High and Ul- 
tra High Energy Cosmic (UHECR) Ray Physics. Recently, 
in particular, the Pierre Auger Collaboration |1 1 and the 
Telescope Array Collaboration |2|, thanks to the excellent 
performance of their hybrid detector arrays, are provid- 
ing us new exciting observations of UHECRs. Although 
these recent results have brought a deeper insight in pri- 
mary cosmic ray properties, still they are largely aflTected 
by the poor knowledge of the nuclear interactions in the 
earth's atmosphere. A calibration of the energy scale in the 
10^^ 10^^ eV energy range accessible to LHC provides 
crucial input for a better interpretation of primary cosmic 
ray properties, in the region between the "knee" and the 
GZK cut-off. 



2 The LHCf experiment 

The LHCf experiment is composed by two independent 
position sensitive electromagnetic calorimeters, located on 
both side of the ATLAS experiment, 140 m away from the 
LHC-IPl interaction point, inside the zero-degree neutral 
absorber (Target Neutral Absorber, TAN). Charged parti- 
cles from the IP are swept away by the inner beam separa- 
tion dipole before reaching the TAN, so that only photons 
mainly from tt^ decays, neutrons and neutral kaons reach 
the LHCf calorimeters. 



Each calorimeter (ARMl and ARM2) has a double tower 
structure, with the smaller tower located at zero degree col- 
lision angle, approximately covering the region with pseudo- 
rapidity 7/ > 10 and the larger one, approximately cover- 
ing the region with 8.4 < r] < 10. Four X-Y layers of 
position sensitive detectors (scintillating fibers in ARMl, 
silicon micro-strip detectors in ARM2) provide measure- 
ments of the transverse profile of the showers. The two 
tower structure allows to reconstruct the jfi decaying in 
two ys, hitting separately the two towers, hence providing 
a very precise absolute energy calibration of the detectors. 
In the range E> 100 GeV, the LHCf detectors have energy 
and position resolutions for electromagnetic showers better 
than 5% and 200yum, respectively. A detailed description of 
the LHCf experimental set-up and of the expected physics 
performances can be found in Ref . |3J . 



3 The single photon energy spectrum 

The LHCf Collaboration has recently published the mea- 
surement of the single photon energy spectrum at 7 TeV p- 
p collisions in two pseudo-rapidity bins |4|. For this anal- 
ysis only a small subset of the 7 TeV p-p collision data 
has been used, corresponding to an integrated luminosity 
of 0.68 nb"^ and 0.52 nb"^ for the ARMl and ARM2 de- 
tectors, respectively. The analysed data have been chosen 
in a particularly clean and low luminosity fill, to minimize 
backgrounds, hence reducing the systematics of the mea- 
surement. 

The main steps of the analysis work-flow with spe- 
cial emphasis to possible implication for the calibration 
of Monte Carlo models used in HECR Physics are sum- 
marised in the following. 

The energy of photons is reconstructed from the signal 
released by the shower particles in the scintillators, after 
applying corrections for the non-uniformity of light col- 
lection and for particles leaking in and out of the edges 
of the calorimeter towers. In order to correct for these last 
two efifects, which are rather important due to the limited 
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Fig. 1. Single photon energy spectra measured by LHCf (black 
dots) for 7] > 10.94 (left) and 8.81 < 77 < 8.99 (right) com- 
pared to the predictions of DPMJET III 3.04 (red), QGSJET 11-03 
(blue), SIBYLL 2.1 (green), EPOS 1.99 (magenta) and PYTHIA 
8.145 (yellow). Top panels show the spectra and bottom panels 
show the ratio of MC results to experimental data. Error bars and 
gray shaded areas in each plot indicate the statistical and the sys- 
tematic errors, respectively. The blue shaded area indicates the 
statistical error of the MC data. Figure from Ref. 14 J . 



transverse size of both the calorimetric towers, the trans- 
verse impact position of showers provided by the position 
sensitive detectors is used. 

Event produced by neutral hadrons are rejected using 
the information about the longitudinal development of the 
showers, which is different for electromagnetic and hadronic 
particles. In addition, thanks to the information provided 
by the position sensitive detectors, events with more than 
one shower inside the same tower (multi-hit) are rejected. 
In order to combine the spectra measured by ARMl and 
ARM2, which have different geometrical configurations, 
in this analysis only events detected in a common pseudo- 
rapidity and azimuthal range are selected: t] > 10.94 and 
A(p = 360° for the small towers and 8.81 < 7/ < 8.99 
and A(p = 20° for the large towers. Figure [T] shows the 
single y spectrum measured by LHCf in the two pseudo- 
rapidity regions compared with results predicted by MC 
simulations using different models: DPMJET III-3.04 0, 
QGSJET 11-03 |6||, SIBYLL 2.1 |7|, EPOS 1.9 and 
PYTHIA 8.145 121. Statical errors and systematic uncer- 
tainties are also plotted. A careful study of systematic un- 
certainties has been done and conservative estimates have 
been taken into account. Further details can be found in 
Ref. |4|. As can be seen from Fig. [T] a clear discrepancy 
between the experimental results and the predictions of the 
models in the whole energy region is present. 



4 Impact of LHCf results on HECR Physics 

The first LHCf results on the photon energy spectra have 
raised attention in the HECR community. As reported in 
the previous paragraph, none of the models agree in the 
whole energy range with the data. Tuning of the models 
are hence needed to describe the Physics of hadronic inter- 
actions at the TeV scale. In order to better understand the 
implication of this measurement for the HECR Physics, 
a collaboration with several MC developers and theoreti- 
cians has started. 



10000 
1000 
100 
10 
1 

0.1 
0.01 
0.001 



I I I I I I I 



DPMJET3 original 
Artificial modification 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



TT° spectrum at E^^^ = lO^^eV 

Longitudinal AS development 



7e+07 
6e+07 
5e+07 
4e+07 
3e+07 
2e+07 
1e+07 



<X^ax>=718g/cm^ 
: <X^ax>=689 g/cm^ 



100 200 300 400 500 600 700 800 900 1000 

Vertical Depth (g/cm^) 

Fig. 2. spectra at E^^ = 10^^ eV for original DPMJET III- 
3.04 model (red) and artificially modified models (gree, blue, ma- 
genta), see text (top). Average longitudinal development of the 
atmospheric shower in the DPMJET III-3.04 model (red points) 
and in the artificially modified model (green points) (bottom). 



As an example, we have artificially modified the DP- 
MJET III 3.04 model to produce a tt^ spectrum that diflTers 
from the original one by an amount approximately equal to 
the difference expected between the different models. Fig- 
ure[2](top) shows the tt^ spectra at Eiab = 10^^ eV predicted 
by the original DPMJET III 3.04 model and the artificially 
modified ones, while on the bottom panel the average lon- 
gitudinal development of the atmospheric shower in the 
original model (red points) and in the artificially modified 
model (green points) is shown. A difference in the posi- 
tion of the shower maximum of the order of 30 g/cm^ is 
observed. 

Fig. [3] shows, as an example, the most recent results 
by Auger Collaboration |10| for the distribution of the 



> variable as function of the energy, which is 
the most used method to infer cosmic rays composition, 
compared with the model predictions for a proton-like (red 
lines) and an Iron-like (blue lines) cosmic ray components, 
respectively. The difference in the < X^ax > distributions 
for the two cases is of the order of 100 g/cm^ hence a 30 
g/cm^ shift is a sizable difference which may significantly 
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Fig. 3. < Xmax > distribution as measured by Auger 1 10 1 (black 
points) compared with the model expectations for a light (red) or 
heavy (blue) cosmic ray composition. The yellow arrow corre- 
spond to the 30 g/cm^ shift, obtained in Fig. 2. 



reflect in the interpretation of HECR data. 

The importance of a direct measurement of the y and tt^ 

spectra by LHCf is clear. 



5 Perspective for a Proton/Lead run at LHC 

The data gathered by LHCf both at 900 GeV and 7 TeV 
center of mass energy are extremely useful for the cali- 
bration of hadronic interaction models used for the study 
of the development of atmospheric showers. Anyway the 
study of the very forward interactions in the pp system 
does not allow to have a complete picture of hadronic in- 
teractions in atmosphere, where mainly nitrogen and oxy- 
gen nuclei are involved, since the eff'ect of the medium can 
not be neglected if precise estimations are requested. The 
LHCf original physics performances can hence be nicely 
improved by including measurements of the particles very 
forward produced in the proton-Nuclei (pA) collisions. In 
principle the ideal way would be to perform the measure- 
ment of proton interactions with nitrogen or oxygen nu- 
clei. This is not possible at the moment neither at the LHC 
nor at other facilities, at least for TeV energies. The LHC 
though, has instead the unique possibility to study the nu- 
clear eff'ects in the forward region of the collisions between 
protons and Lead nuclei at the TeV scale, corresponding to 
a proton energy in the LAB frame exceeding 10^^ eV. Pre- 
vious measurements done at other facilities ifTTIl reported a 
marked reduction of the production cross section values for 
low transverse momentum events in the forward region of 
ion interactions, with respect to what could be expected by 
considering the events as binary collisions between nucle- 
ons. These results, normally reported as the ratio between 
cross sections measured in p-ion (or ion-ion) and in pp col- 
lisions (that is the Nuclear Modification Factor, NMF), are 
currently explained with nucleon screening eff'ects and sat- 
uration of partonic densities, but so far we have only a poor 
knowledge of the mechanisms involved, especially at very 
high rapidity where no measurements currently exist. 

To experimentally investigate these important physics 
aspects, the LHCf collaboration has proposed to the LHCC 
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Fig. 4. Multiplicities of photons and neutron hits on the proton- 
remnant side expected from the EPICS and DPMJET-III models 
for p-Pb collisions at 4.4 TeV center of mass energy in the proton 
remnant side. 



to take data during the next proton-Lead run at LHC, cur- 
rently foreseen for the end of 2012. We have submitted to 
LHCC a Letter of Intent 1 1 2 1 , detailing the physics case for 
this type of measurement, and showing that the existing de- 
tector configuration (Arm2 in particular) is well suited for 
this purpose. Some interesting results are shown in the last 
part of this section, based on the simulations performed 
using EPICS and DPMJET-III hadronic interaction mod- 
els to describe the proton-Lead collisions. A total of 10^ 
events have been generated for each model, considering a 
3.5 TeV proton hitting a 1.38 TeV/nucleon Lead ion with a 
null crossing angle and random impact parameter. The to- 
tal energy in the nucleon-nucleon center of mass frame is 
approximately 4.4 TeV. The most significant results are ob- 
tained by measuring particles produced in the proton rem- 
nant side, that is more simply accessible from the experi- 
mental point of view; in this region smaller multiplicity of 
neutral particles are expected, as can be seen from Eig. [4j 
that reports the expected number of photons and neutrons 
hitting both the Arm2 towers. 

Approximately 1% of the events have one single hit, 
and less than a fraction of 10"^ of the events have two hits 
on a single tower. Therefore we do not expect particular 
problems both during data taking and analysis phases. Ad- 
ditionally, we want to stress that the absolute energy scale 
can be directly checked also in the p-Pb run with the data 
by the direct tt^ reconstruction method, by separately iden- 
tifying the two photons of the tt^ decay in the two separate 
towers, as has already been demonstrated with the 7 TeV 
pp data |[T3]| . 

Eigure[5]shows the expected photon spectra in the small 
(top) and large (bottom) tower, as obtained with the EPICS 
and DPMJET-III models, for 10^ inelastic p-Pb collisions. 
A clear and significant diff'erence between the two mod- 
els can be seen from this figure. In particular we find that 
the DPMJET-III model has a harder spectra than the EPOS 
model. The diff'erence in the slopes is evident and LHCf 
will easily discriminate between the two models. 



6 Future activities 

New analyses of data collected in 2009-2010 runs are in 
progress, in particular the measurement of the single pho- 
ton spectra at 900 GeV is almost finalised as well as the 
reconstruction of tt^ spectra. 
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Fig. 5. Expected energy spectra of photons hitting the Arm2 de- 
tector on the proton-remnant side small tower (top) and big tower 
(bottom). The statistics shown in the plots correspond to a data 
taking period of few hours at the nominal 10^^cm~^s~^ luminos- 
ity. 





Fig. 6. LHCf Pt vs energy acceptance region covered by type I 
(left) and type II (rigth) tt^ reconstruction 



The possibility to reconstruct also tt^ from pair of pho- 
tons impinging in the same tower (type II) in addition to 
the TT^ reconstruction from pair photons impinging one in 
each tower (type I) is under investigation. As can be seen 
from Fig. |6] the reconstruction of type II tt^ events gives 
the possibility to explore a different pr region: while type 

I events are mainly characterised by pair of photons with 
large opening angle and dominates at low energies, type 

II TT^ events are characterized by tight photon opening an- 
gles and dominates at higher energies. The combination of 
type I and type II reconstruction thus enlarges our physics 
acceptance. 

In the meantime, the LHCf Collaboration is also work- 
ing on the upgrade of the detector to improve the radia- 
tion resistance in view of the 14 TeV p-p run, currently 
foreseen in 2014. The scintillating part of the detector will 
be replaced with GSO slabs. Beam test results |[T4l have 
demonstrated that the upgrade of the calorimeter with GSO 
will enable LHCf to sustain the radiation level foreseen 
in the 14 TeV run. Additional improvements in the front- 
end electronics of the silicon position sensitive layers of 
ARM2 detectors as well as an optimization of the layout 



to improve the stand-alone silicon energy resolution are on 
going. 



7 Summary 

LHCf has measured for the first time the single photon 
energy spectra in the very forward region of 7 TeV p-p 
collisions at LHC. None of the hadron interaction models 
mostly used in HECR Physics show a perfect agreement 
within the statistical and systematic errors with the data. 
New analyses and additional studies are in progress and 
the detector upgrade for the p-p run at 14 TeV is on going. 
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